Superradiance in spin-J particles: Effects of multiple levels 
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We study the superradiance dynamics in a dense system of atoms each of which can be generally 
a spin-j particle with j an arbitrary half-integer. We generalize Dicke's superradiance point of view 
to multiple-level systems, and compare the results based on a novel approach we have developed 
in [Yelin et al, arXiv:quant-ph/0509184 . Using this formalism we derive an effective two-body 
description that shows cooperative and collective effects for spin-j particles, taking into account the 
coherence of transitions between different atomic levels. We find that the superradiance, which is 
well-known as a many-body phenomenon, can also be modified by multiple level effects. We also 
discuss the feasibility and propose that our approach can be applied to polar molecules, for their 
vibrational states have multi-level structure which is partially harmonic. 

PACS numbers: 42.50.Nn, 42.50.Ar, 33.20.Tp 



I. INTRODUCTION 

Quantum many-body physics has been one of the most 
attractive areas for decades along with the remarkable 
advents in the fields of ultracold atomic and molecu- 
lar systems and quantum optics. These systems not 
only provide an excellent testbed to study quantum na- 
ture of various many-body phenomena such as Bose- 
Einstein condensation, superfluidity 1-31, quantum mag- 
netism, and quantum phase transitions 4j-[7| , but also in- 
spire the implementation of quantum machinery such as 
quantum simulation (sl-fisj and quantum computing 
|l6|. In many ways, quantum many-body effects are "ex- 
otic" compared to their classical counterparts, and even 
to quantum single-body physics, mainly due to particle 
statistics and indistinguishability of particles. The cir- 
cumstances can become even more complicated when an 
ensemble of particles interact cooperatively which results 
in higher-order nonlinear effects. Superradiance, usually 
representing an iV^ enhancement of the radiation inten- 
sity due to coherent decay of a dense sample consisting of 
N excited atoms, is one important example which can be 
understood qualitatively through particle indistinguisha- 
bility and symmetry arguments without the need for con- 
sidering particle statistics. This phenomenon was first 
predicted in 1954 by R. H. Dicke [l3|, who pointed out 
that the radiative properties of an excited atom can be 
very different just because other atoms are present or not, 
given that their distance is much smaller than the wave- 
length of the radiation field even if the particle wavepack- 
ets do not overlap and no direct interaction is present. 
Since then such cooperative effects have been intensively 
investigated both theoretically and experimentally |18l - 
H |. Recently, superradiance re-catches one's attention 
in the context of a Bose-Einstein condensate coupled to 
a cavity 25, 26], alkaline-earth-metal atoms ^7^, Ryd- 
berg atoms [2J, [51] , and quantum dots [1^ , as well as its 
strong connections to quantum information through the 
so-called Dicke states. Such states are fully symmetric 
states by particle permutation, and mostly serve as the 
main stage during the superradiance process. 



Traditionally, superradiance deals with two-level 
atoms or other spin-Y2 systems, being first excited, that 
decay cooperatively. It is natural as a next step to 
consider particles with larger spins, i.e. systems with 
a multiple level structure. Examples include the near- 
harmonic level structure of low-lying vibrational states 
in molecules. Generally speaking, multi-level structure 
brings up more complications to the radiating system. 
For example, an excited atom in a higher level is still 
excited after emitting a photon; the atoms and photons 
from different-level transitions can further cooperate and 
modify the overall emission behavior. In order to study 
multi-level effects, we re-consider Dicke's point of view of 
superradiance as the starting point by first assuming the 
system to be point-like and fully symmetric. 

However, Dicke's picture is only qualitatively correct 
and insufficient to describe real situations, where the ac- 
tual arrangement of particles, the sample's finite size, 
and dipole-dipole interactions play a role. Microscopi- 
cally, single atoms build up inter-atomic coherence due to 
virtual photon exchange caused by dipole-dipole interac- 
tion and form many-body states such as the Dicke states 
[Tsl [l9| . The coherence can be breached when the ge- 
ometry of particle arrangement introduces inhomogene- 
ity such as dipole-dipole interaction between each pair of 
particles. This leads to dephasing effects and therefore 
Dicke's picture fails to be valid. To characterize how the 
"finite size" influences superradiant behavior, a param- 
eter, cooperativity C ^ MX^ is introduced, with M the 
number density and A the wavelength of the transition 
field. One then expects that superradiance is observable 
for C 3> 1 and be suppressed for small C. Our previous 
study [28i] has further suggested a more accurate esti- 
mation that the criterion of observing superradiance is 
approximately given by AfX^d (d is the sample size), in 
agreement with [24], where such a factor is found to be 
an essential one that determines whether cooperative ef- 
fects dominate or not. To take into account the realistic 
arrangement of our particle systems, we use a novel for- 
malism that considers only two probe particles, treating 
the spread of environment atoms in the mean-field ap- 
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Figure 1: (color online). The energy level structure for A'^ = 3 



tion plays no role as well as the fermionic or bosonic 
nature of the atoms. Suppose that the transitions are 
induced by dipoles through the interaction Hamiltonian 
V = — ^j^Pi ■ E{fi]^, where pi is the dipole operator of 
the ith atom and E{'ri) is the local field at the coordi- 
nate fi. Under the long- wavelength assumption of the 
field for a given small system size, the spatial depen- 
dence can be eliminated. Therefore V = —E ■ J^iPi — 

~J2^l=x,y,zP^^i^^^^^ + ^■^■) rotating wave ap- 

proximation. Here is the dipole moment magnitude 



three- level atoms {j = 1). There are {2j + if = 27 levels, of an atom in the fi direction, is the positive (nega- 



including the fully symmetric 2Nj + 1 = 7 states (the Dicke 
states) . 



proximation and then take an average over all possible 
particle pairs [13, [Mj . This approach enables us to write 
down an effective master equation, retaining the degrees 
of freedom of two-body coherence which can be regarded 
as a projection of the many-body coherence in the orig- 
inal system. This method has been proven to show a 
good agreement with ongoing experiments with Rydberg 
atoms [28| . In this paper we further apply this formalism 
to spin- J systems. 

This article is organized as follows: In Sec. [n]we dis- 
cuss the original picture proposed by Dicke and gener- 
alize the idea to multi-level srotems. Sec. IIIII sketches 
the formalism developed in [sO, HH and summarize the 
governing equations. We then apply this method to mul- 
tiple levels and show the results in Sec. IIVI There we 
also discuss the differences from the Dicke model and in- 
vestigate the significance of many-body and multi-level 
correlations. In Sec. |V]we further consider the thermal 
Doppler broadening and calculate the marginal condi- 
tions that superradiance can tolerate. Finally, in Sec. IVII 
a dipolar molecular gas is discussed as an example, for 
which we consider the vibrational states and investigate 
the superradiance effects from its vibrational states. 



II. DICKE SUPERRADIANCE 

To gain qualitative understanding of the Dicke superra- 
diance, we start with considering an ensemble of (2j -1-1)- 
level atoms confined within a small region with its size 
much smaller than the wavelength of the radiation field. 
In this limit, the particles are indistinguishable viewed 
by the field and must be regarded as a whole quantum 
object. To emphasize that the collective radiative be- 
havior is governed solely by many-body effects, we do 
not assume any instantaneous, i.e. non-radiative, in- 
teraction between atoms. The inter-particle spacing is 
large enough so that the overlap of particle wavepack- 
ets is negligible. In other words, the exchange interac- 



tive) frequency component of the field, and = af^ 

with = X]m=-,7 + 1| and (t+ = {a^y . Note 

that V does not break the permutational symmetry of 
the particles. If all the atoms are initially excited, time 
evolution will only take the state of the system around 
the fully symmetric manifold, whose eigenstates are usu- 
ally called the Dicke states (See Fig. [TJ [13, [H): 



|J,Af) 



' (J + M)! 
(2J)!(J-M)! 



{j-y-''\j,j), 



(1) 



where J = Nj is the total spin of N spin-j atoms and the 
integer M denoting the level index can only go from J 
through — J; the total spin ladder operators J"^ — J^iJi" 
satisfy J±|J,A/) = y/ J{J + I) - M{M ± 1)\J, M ± 1) 
with each satisfying an analogous relation within the 
ith atom. The emission rate is then given by = 
J2m PmWj{M), where pM is the probability for the state 
being at the Mth level, and the associated collective de- 
cay rate is Wj{M) = ^{D^ D~) jm with 7 denoting the 
bare rate in free space. 

For spin-i/2 particles, the ladder operator happens 
to be, up to a constant factor, equivalent to the dipole 
operator, = |5)i(e| and af — \e)i{g\. This connection 
makes it straightforward to obtain {D~^D~)jm = {J + 
M){J - M + 1) [li|. For spin-j > 1/2 atoms, the spin 
and dipole operators are no longer parallel. To obtain 
an explicit relation in this case, we take the mean-field 
assumption and get 
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IJM 



= N{a+ai).JM + N{N - l){a+a^).,Km 

where {(J^<Ti)jm and {<J^(J2)jm can be further ex- 
pressed in terms of the Clebsch-Gordan coefficients: 



(a+a^) = 1 - {j, (N ~ l)j; ~j, M + j\J, M)' (3) 



and 



3 



CTl (Jo = 



J2 [O'. - mi, M ~ mil J, Af) (j, (A^ - 2)j; ma, M - mi - maj J, M - mi) x 

mi ,7712 

(j, (iV - l)j; mi - 1, M - mi + 1| J, M) (j, (iV - 2)j; ma + 1, Af - mi - maj J, Af - mi + 1) . (4) 
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Figure 2: (color online), (a) (l)+l)^)./Af as a function of Af — 
J and (b) the emission intensity per particle Jem of different 
spin j's for A'^ = 10 atoms; (c) Overall emission curves for 
different spin j's when J — Nj — 15 is fixed. Note that Jem 
and W are in units of 7^010 with the energy level spacing hcjo- 



The equation of motion now reads 
Pm=j = -Wj{J)pj, 

PM<J = -Wj{M)pM + Wj{M + l)pM+l. 



(5) 



The emission curves of different j are shown in Fig. 
[2Jb), for which we calculate the intensity per particle 
4^ = W{t)/N with iV = 10 by evolving Eq. ©. One 
can observe that every curve shows a different degree 
of superradiance behavior, i.e., the intensity grows and 
maximizes in a short period of time. As j increases, the 



peak intensity becomes higher. This implies that the 
radiation enhancement not only comes from many-body 
effects but also from multiple levels. This can be also seen 
in Fig. [2{a), where the enhancement factor {D^D^)jm 
as a function of Af — J is plotted, consistently explain- 
ing the higher emission rate for larger j. One feature 
worth noting is that all j-curves in Fig. [HJb) share the 
same "growing" behavior. This is also suggested by the 
{I)^I}~)jM curves in Fig. HJa): For we start the ra- 
diation process from the fully excited state (the highest 
collective level), the dynamics is dominantly determined 
by the population flows associated with a few higher lev- 
els. The rates are proportional to the enhancement fac- 
tor {D~^D~)jM- For various j, we find in Fig. ^a.) that 
these curves coincide on the left hand side, corresponding 
to |Af — J| small (highest levels). Another noticeable fea- 
ture is that the (f)+l)~) jm curves develop a plateau as j 
increases and the highest value is found to be bounded in 
the large-j limit. It can then be expected that the peak 
intensity value should also have a bounded value for very 
large j. We will find this observation still true when we 
use a more sophisticated approach. More details will be 
discussed in Sec. IIIIl 

Finally, we compare the overall radiation by differ- 
ent spin-j atoms, keeping the total spin J = Nj fixed. 
Among these cases, they share the same J-Bloch sphere 
and therefore might be expected to have similar behav- 
iors. However this is not true, as we can see from Fig. 
^c). Smaller j cases have faster and more intense burst 
of emission while larger j have smoother emission rates 
and longer tails. This is because {D^D~)jm, as deter- 
mined by Eqs. ^ and is also dependent on A^, not 
merely dependent on the total spin J. 



III. EFFECTIVE TWO-BODY FORMALISM 

We would like to emphasize that dipole-dipole inter- 
action plays a crucial role and is responsible for both 
real and virtual photon exchange. As a consequence, the 
system builds up inter-particle coherence while decaying. 
As dipole-dipole interaction is contained in Dicke's pic- 
ture in a sense that the spin flip-flops count, this picture 
treats the whole system as a point-like object so that the 
inter-dipole coupling is considered uniform. In an ac- 
tual laboratory setup, Dicke's picture is usually an over- 
simplifled view because a real sample always occupies 
a finite size and sees a finite wavelength of the radia- 
tive field. The spatial arrangement of particles usually 
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breaks permutation symmetry. (Or more precisely, each 
particle sees different dipole-dipole couplings to all oth- 
ers.) The nonuniform coupling leads to dephasing of the 
Dicke states, resulting in population leakage out of the 
fully symmetric manifold. Furthermore, dipole-dipole in- 
teraction also causes other effects, e.g. frequency chirp- 
ing, for which each Dicke state | JM) can be dipole-dipole 
shifted differently so that the emission frequency becomes 
variable over time [loj . Superradiant behavior becomes 
more complex (and less pronounced) when these effects 
are not excluded. In order to better describe practical 
situations, we need to go over the microscopic details 
of atom- field interactions. The calculation, however, be- 
comes intractable when the number of particles increases 
typically for > 10. In f2l, [sO, |3l| , we circumvent this 
difficulty by explicitly writing down the master equation 
of motion for only two probe atoms, taking average over 
the background atoms and tracing out the field variables. 
We also assume that the field instantaneously interacts 
with the whole ensemble, ignoring the retarding effects 
due to the finite size. (This can be justified because the 
characteristic time of propagation d/c (~ 10~^^ sec for a 
sample of size d ^ 1 mm) is usually much smaller than 
any other decay timescales.) We summarize here the 
main results and leave the details of the derivation to 
the Appendix. The relevant two-body master equation 
is given by 

r,. 



r 



l^[[pa],a-] + [a]^a-p\), (6) 



- E 

where p is the two-body density matrix with dimension 



(2j + If X {2j + 1)2, 7 



is the free-space spon- 



taneous decay rate, T = Ta is the single-particle induced 
pump/decay rate, and f = [i ^ j) denotes the two- 
particle damping rate responsible for the atom-atom cor- 
relation. The mean-field approximation with the second 
order correction in fields gives the self-consistent form for 
the induced rate: 
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where p^^^ = \ 2^^»M denotes the reduced single- 

particle density matrix and Pab;cd = jii^^j c||p|6, c?) -I- 
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Figure 3: (color online). The temporal profiles of (a) the 
emission intensity per particle (in units of ^huo) and (b) the 
induced single-atom pump/decay rate F and two-atom cor- 
relation damping rate F (in units of 7). In all cases we use 
C = 10 and £1 = 10. 



{c,a\p\d,b)]. The factor i comes from averaging the in- 
terchanging of two particles. Note that interchange sym- 
metry requires pab.cd = Pcd,ab and /o*^.^^ = pba-dc- The 
cooperativity parameter is defined as C = /u'^] 
Q = ujd/{2c) characterizes the system size d in terms of 
the radiation wavelength. These results are based on as- 
sumptions that there is no external field and hence the 
generated field has to be on-resonant with the transition 
frequency. Function /(C, p) = [(((■— l)e''-^'"^s n") 4-(npi- 



singy]/ie + 6'? 



-cosg) +[ge 
i~2_^g2 for large C, and g. If no ther- 
mal broadening is assumed, we have ( = ^Cg ^_^_^^^ V{t). 
When the Doppler broadening needs to be considered, 
the fields allow to be detuned and these quantities must 
be averaged. More details will be discussed in Sec. |Vl 



IV. RESULTS 



A. Emission and decay rates 

By Eq. (O we are able to solve for the temporal emis- 
sion rate curve. Fig. [3{a) shows the radiation intensity 
per particle Jem with C = 10 and g = 10 for different 
spin-j species, where /cm = fu^o I]m=-j + 
Here we take Pjjjj = 1 and set for all other density 
matrix elements as the initial state. It can be seen that 
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for each j the radiation intensity reaches a peak, giving 
a strong evidence of superradiance. Those curves fol- 
low roughly the same intensity profile in the beginning. 
The maximal value of intensity first grows as j increases 
from 1/2, and then stops growing when j > 2. The time 
of reaching the peak intensity also converges to a fixed 
constant imax in the large j limit. This also has been 
observed in the Dicke superradiance picture. Here we 
give an intuitive explanation as follows: When the sys- 
tem starts to relax from the state with all atoms initially 
excited to the highest level, only a few highest levels are 
involved in determining the radiative behavior during the 
early stage. Even for a very large spin particle, who has 
a huge multi-level structure, those levels lower than the 
first few have not been populated yet and hence do not 
have contributions. The time evolution of the decay rates 
is also plotted in Fig. [3)Jb). Note that at t = the di- 
agonal decay rate r(0) determines the initial emission 
intensity, followed by a sharp growth and hence resulting 
in intensity peaks. In the mean time, the off-diagonal 
r emerges and mixes single-body states. Different from 
the Dicke picture, where we choose the eigenbasis to be 
the symmetric states constructed by a giant spin object 
J = Nj, here we use products of single particle states as 
the eigenbasis, allowing the degrees of freedom of popu- 
lation being transferred to asymmetric levels. Note that 
the dipole-dipole interaction is built-in in our formalism 
and is responsible for these effects. Consequently, the 
superradiance enhancement with j, when characterized 
by the growth of peak intensity, in more realistic cases 
cannot be as large as predicted by the Dicke model. On 
the other hand, since the asymmetric levels have lower 
or vanishing decay rates, the occupation of these levels 
modifies the tails of the emission curves. In some cir- 
cumstances, the energy is trapped. Such effects cannot 
be described by the Dicke model. 

A few remarks are placed here regarding the connec- 
tion of cooperativity and superradiant curves. As we in- 
crease C (g) while fixing g (C), the emission peak intensity 
per particle increases proportionally while the time scale 
of the initial intensity burst is inversely proportional to C 
(g). This is due to the "many-body enhancement", as we 
have discussed using Dicke's picture. Such features are 
commonly observed even in the original two-level sys- 
tems. Furthermore, the emission curves are found to 
be similar when Cg is kept the same (not shown). This 
can also be seen analytically (see Appendix). This sug- 
gests that Cg ^ AfX^d be the relevant factor that deter- 
mines the primary superradiant behavior while the de- 
tailed emission curves, however, still slightly depend on 
C and g separately [ll,!!!]. 



B. Significance of atom-atom coherence 

Cooperation of many-body states is crucial to superra- 
diance. The beauty of this formalism is that we retain the 
accessibility to atom-atom correlations under the frame- 
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Figure 4: (color online). Evolution curves for emission inten- 
sity when the off-diagonal terms are fully considered (full), 
partially removed (with both (i) and (ii), or only with (i); see 
text), or entirely removed (no off-diag) for (a) spin-i/2, (b) 
spin-1, and (c) spin-9/2 particles. Other parameters are the 
same as Fig. O Inset of (a): We plot three curves of su- 
perradiance (SR, blue solid), amplified spontaneous emission 
(ASE, red dashed), and single-particle free-space spontaneous 
emission (black dotted) for comparison. 



work of the mean field approximation. Here we inves- 
tigate the role of many-body correlations, which, in our 
case, are contained in the off-diagonal terms of the two- 
body density matrix. To see this, this method allows us 
to manipulate these off-diagonal terms and evolve Eq. 
([5]). The results will then be compared. Note that these 
off-diagonal terms have the form Pa,a+m;b,b-m- In the 
spin-i/2 case, the only possibility is Peg,ge- Fig- Hla) 
shows the emission curve for C = 10, g = 10, and 
Peg,ge = at all time. This curve is now found to be 
monotonically descending, signaling mere amplified spon- 
taneous emission (ASE) instead of superradiance (SR) 
because of apparent lack of an intensity peak. Such 
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monotonicity is shared by larger j cases (Figs. IHJb) and 
(c)) when aU ofT-diagonal terms are set to zero. The 
reason is clear: Without atom-atom correlation, the den- 
sity matrix is reduced to a single-particle description and 
therefore no cooperative effects are observed. On the 
other hand, for larger j atoms the off-diagonal terms does 
not only concern atom-atom correlations from the same- 
level transitions but also involve that from different-level 
density matrix elements. In the following we try to dis- 
tinguish the importance from three kinds of coherence 
terms: (i) the same-level coherence Pa.a+i:a+i,m (h) the 
cross-coherence Pa,a+i-b+i,b for a ^h, and (iii) the higher- 
order coherence Pa,a+m:b,b-m for m > 2. For example, 
in Figs, m^b) for spin-1 and (c) for spin-9/2 atoms, we 
plot the emission curves corresponding to all off-diagonal 
terms being dropped out, and (i), (ii), and then (iii) being 
added back to the system. It can be found that by inclu- 
sion of (i) and (ii) the system has already behaved like 
the actual dynamics, indicating that the higher-order co- 
herence is negligible in determining the evolution of the 
system. However, if only (i) is included, although the 
intensity enhancement can still be observed, the details 
of emission profile have discrepancy to the actual behav- 
iors. The distinction becomes even more obvious when 
j gets large, as we see in Fig. HJc). The fact that the 
cross-coherence terms must be taken into consideration 
implies the interferences due to "cross-level" transitions, 
i.e., the transitions of the same energy difference, but 
not from two definite levels, have some kind of "multi- 
level" contributions to the superradiance, in analog to 
the many-body effects. 



(I43I44I42P and Green's functions (gS]) need to be deter- 
mined in an averaged manner because they contain the 
effects from all slivers. For convenience we use a nota- 
tion that Q — Q{5)fD{5)d5 to represent an Doppler 
averaged quantity. Consequently, 
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where U{z) = /z°° is the scaled comple- 

mentary error function and zo = ^^y^A^'*^ ' '^o ~ 
^V{t)^y^U{izo) while Yo ^ lo ^ uj/c keeps the 
same. Finally we have 
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_2J_c2//(C,e)Re[[/(*zo)]r(i) (16) 



f (A) = -^/(C, g)Re[U{izo)] x 
[3CgA{t) + 2C2/r(i)] , 



(17) 



where 



V. DOPPLER BROADENING 

When a hot gas is considered, the energy difference 
seen by moving atoms varies. In this section we consider 
the loss of coherence due to the Doppler effects. Suppose 
that the thermal gas is described by a Gaussian distribu- 
tion function: 
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(18) 
(19) 



Each sliver, depending on its location at the thermal dis- 
tribution function, sees these quantities associated with a 
specific Fourier component A. The overall decay is then 
equivalent to averaging F and F according to the thermal 
distribution function, i.e., 



where S is the Doppler shift for an atom and Ajj is the 
characteristic width of this distribution. In order to take 
this average into account we need to go to the deriva- 
tion summarized in the Appendix. Note that the fre- 
quency difference between the field and the Fourier com- 
ponent A in Eqs. I|40l4ip is now modified to A — (5. 
Here we also assume that the velocity distribution will 
not be affected by atoms' recoil when photons are emit- 
ted. The idea is then to divide the whole system in the 
frequency space into a series of small slivers according to 
the detunings they see individually. For each sliver, it 
can be expected that the integration as in Eqs. I|4l)l4ip 
should lead to the same spatial dependence. The only 
modification is that the source and retarded functions 
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-e-^'/(2Ai)ry(A). 



(20) 



Eq. (|2Up is solved numerically. We then calculate the 
corresponding emission curves and show them in Fig. [5l 
With the Doppler broadening, it is clear that the super- 
radiance behavior is suppressed as the Doppler width A d 
increases. This is due to frequency mismatch and there- 
fore part of the atoms loses track of coherence and decays 
more independently. But the superradiance peaks are 
still observable within a certain range of A^j, until it be- 
comes too large and kills the peaks. This is understand- 
able because for a small-width distribution, there are 
still sufficient atoms located within a frequency-matching 
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Figure 5: (color online). Doppler-broadened emission inten- 
sity curves for (a) spin-i/2 and (b) spin-9/2 particles. The 
inset figures show, for each case, two curves slightly smaller 
and larger than the marginal Doppler width Am- Here we use 
C = 10 and g — 10 for both cases, and get (a) Am/7 — 433 
and (b) Am/7 = 1650. 
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Figure 6: (color online). The marginal Doppler width Am as 
a function of the cooperativity for (a) spin-i/a and (b) spin- 
9/2 particles. (Here we take q = 10). The cross dots (blue) 
represent our calculated data; the solid lines (red) are best 
fitting power-law curves with the exponent around 2. 



regime. In order to characterize the existence of superra- 
diance when the match and mismatch parts compete, we 
define a marginal Doppler width A„j beyond which the 
peak value no longer surpasses the initial intensity. Obvi- 
ously A„j depends on the cooperativity parameter C, and 
in Figs. El[a) and (b) we plot the relations for spin-i/2 and 
spin-9/2 particles, respectively. We find approximately a 
quadratic dependence A^ oc C^. Note that when the co- 
operativity parameter increases, not only the total num- 
ber of particle increases but the inter-particle spacing 
decreases accordingly. This enhances the dipole-dipole 
interaction (of the order of magnitude p^/(27reor^) oc C) 
which is responsible for the superradiance for an addi- 
tional factor C. As a result we have the tolerance A™ 
with a quadratic dependence on C rather than a linear 
dependence. 

VI. MOLECULAR VIBRATIONAL STATES 

One direct example for multi-level structure is vibra- 
tional modes of polar molecules where the deeply bound 
potential can be well approximated by a harmonic one. 
The number of low lying eigenstates that are quasi- 
equally spaced energy levels can usually be up to a few 
tens. These particles are thus analogous to large "spin" 
particles. Take a typical example of heteronuclear di- 



atomic alkali molecules [s^]: the state X^T,'^ for LiCs 
has an averaged energy spacing wo « 27r x 5 THz. For 
a sample of LiCs molecules with a density A/" ~ 4 x 10^ 
cm~'^, such energy spacing corresponds to cooperativity 
C ^ 20 3> 1. The transitional dipole moment between 
two adjacent vibrational states is about 5 Debye, and 
therefore the single-particle spontaneous emission rate 
7 ~ sec~^. We then expect within this parameter regime 
that superradiance intensity peak can be observed in a 
timescale of mini-sec while Fmax ^ IOOO7. 

The cascade relaxation of excited population from 
higher to lower levels is a reminiscence of motional 
cooling. When the cooperative effect comes into play, 
the down-ladder process will be accelerated because the 
stimulated decay becomes dominant while superradiance 
takes place (without other pumping processes such as 
thermal excitation). As we have pointed out, the in- 
duced rate can be a few orders of magnitude by increas- 
ing the number of particles and hence the cooperativity. 
This suggests a scheme of "superradiance-assisted cool- 
ing". Such scheme may be an alternative to cool vibra- 
tional states of molecules, which, generally speaking, has 
previously been obstructed by the fact that such states 
are only weakly optically coupled. 
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Figure 7: Schwinger-Keldysh contour 



APPENDIX 
A. Two-body master equation 

Although this method of the effective two-b ody de- 
scription has been discussed in great details in [sO, HH , 
we summarize in this section the derivation for the for- 
malism for completeness. We start with the microscopic 
Hamiltonian that reads 



Ho 



-Yp,-{£{f,,t)+E{r,,t)). 



(21) 



i=l 



Here, we separate the field into two parts: the exter- 
nal classical driving field £{ri,t) and the induced local 
field Eiji, t). When a uniform dense gas of atoms is con- 
sidered, it is reasonable to assume that every atom sees 
the same field and the same background due to other 
atoms. On the other hand, in order to take into account 
atom-atom quantum correlation we need to retain ade- 
quate degrees of freedom involving at least two particles. 
Our proposal is therefore to write down an effective de- 
scription for two probe atoms in which all other atoms' 
contribution will be averaged in the mean-field sense and 
appear as parameters in the two-body description. In Eq. 
(I^Tjl V is the interaction of the probe atoms (1 = 1,2) 
with the field and will be treated as a small perturbation. 
Hq consists of the unperturbed atomic and field Hamil- 
tonian as well as the contribution from the background 
atoms. In the interaction picture, the evolution operator 
is then given by 



Si{t) = Tcxp 



Vj{t')dt' 



(22) 



where T is the time-ordering operator. We here introduce 
the positive and negative components x{t) = x~^{t) + 
x~{t)^with a:=^(i) = x(t)e=F''^*, where e 
and x(t) is the slowly varying (compared to the inverse 
of the radiation frequency uj~^) amplitude of the corre- 
sponding quantity. Further note that pf^ = p^(J~ and 

Pin ~ Pfj^'^t ■ -'^^ ^^"^ rotating wave approximation, the in- 
teraction becomes Vi{t) ~ J2ip,{ptiii-^JI +^;7) +h-C.]. Eq. 
can be cast in the framework of Schwinger-Keldysh 



formalism, in which Vi{t) — > V{f) and Si{t) 



Sc = Tc exp 



(23) 



where C denotes the Schwinger-Keldysh contour as 
shown in Fig. [71 and Tc is the contour-oriented time- 
ordering operator, i.e., along the upper branch of the 
contour Tc is the normal time-ordering operator while 
along the lower branch is the inverse time-ordering oper- 
ator. To prevent possible confusion we denote the "time" 
parameter along the Keldysh contour with a check sign. 
We then trace out the degrees of freedom of the fields 
and the background atoms, which leads to an effective 
evolution operator: 



ocff 



(•Sc) field 

exp I M dT^Y^ [p+ {f)£l^ {n , r) 



+ P^^'^)^Un,r) 



(24) 



1 



+ Ptn (^1 )Cv (n , h )p+ (r2 ) j L 



where £^ {vi, f ) = £^ {'ri,T) + {E^ {ri,f)) is the local field 
seen by the probe atom, and the Green's function of the 
interacting field 



AAijiy(n, 7^2) 



{{TcE;in,T,)E+if,,f2))) (25) 
a,,jAruh) = {{TcE+{n,h)E-{f,,h))). (26) 

To get Eq. ([M]! we have used 



(Tcexp[si]) — exp 



m! 



(27) 



with ((•)) denoting the cumulant, which is defined by 

P» = {A) 
{{AB)) = {AB)-{A){B), 

where A and B are operators. We also set the higher- 
order cumulants {{E™)) = for m > 2 by assuming that 
the radiation field is Gaussian. The two-field Green's 
function -Di^ijv('ri, f2), depending on the order of fi and 
T2 on C, has four possible forms: 



D 

D1 



{{TE-in,n+)E+ir,,T2+))) (28) 

{{T-'E-{r,,r,^)E+ir,,T2-))) (29) 

{{E;ir,,n^)E+{r,,T2+))) (30) 

{{E+{r,,T2-)E;^in,n+))). (31) 
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The other Green's function C^^jv has similar relations. 
By the subscript "+" or "-" we denote the upper or lower 
branch for t, respectively. In Eq. those terms like 



/OO pOQ 
dri / dT2X (32) 



where k^ib = 1 for A = B and k^ib = — 1 for A 7^ _B, 
and the subscripts A, B being placed with a~ and cr"'" 
emphasizes that the operators must be in order accord- 
ingly on the Schwinger-Keldysh contour C. We change 
the time variables for the field correlation such that 

T2 = T + t'/2 



tf^jA'T-.T') with Ti = T - t' /2 and 
After some math, we reach 



r - foo 

S^J = T.exp / dr[arir+)£Z^ir,,r)~arir^)£-^in,r)+anr+)£U^^^ 

^^ r.,, (c.,r) _ 2ar(r_)a+(r+)] 

, (w, r) + (w, r) r ^ _ ^ _ 0+/' ^l 

dr—i^ [aJ{T+)cr.^ (t+) + (Tj(r_)cr,, (t_) - 2(Tj(r_)o-^ (t+)J 

DO 

dTHi^^j^{u},T) [crr(r+)cr+(T+) - a-r(r_)cr+(r_) - cr+(T+)CT7 (r+) + (T+ (T_)cr," (r_ )] 

dTh,f,,j^{uj,T)[cr+ {T+)a- {t+) + a+ {T_)a^ {t_)] L 
1 J 

where we have introduced these quantities: 



i 



-00 
00 



^^^I..,jv{T, w) 



dr'{{E-{n,T)E+{r,,r + r')}}e^-^\ 
dT'{[E+{r,,T),E-{r„T + T')])e'-^', 



/ dr'{{{E-{n,r)E+{r„r^r')))e-^-^' - ((i?; (f, , r)£;+(f„ r + r')))e-' }, 



(33) 

(34) 
(35) 
(36) 



2h 



dr' { { [E- in , r) , E+ (r, , r - r')] ) ' - ( [i?; (f, , r) , E+ [r, , r + r')] ) e^^' } . (37) 

I 



From Eq. p3p we extract the effective master equation: as the induced pump and decay rate and term 7 as the 

spontaneous decay rate inside the atomic medium. Term 



z— 1,2 /^i,^ 

r. 



E E^(K^,-]] + K,-1p]) 



- E E '^"'^^r"" (H-^n + K.<p])- 

2, J — 1,2 

Note that p is the two-body density operator so is a 
(2j -I- 1)2 X (2j + 1)2 matrix. Now we identify term T 



h as shown in Eq. (|37p corresponds to the Lamb shifts, 
which are somewhat irrelevant for our current considera- 
tion and is therefore absorbed to the unperturbed Hamil- 
tonian Hq; term H as shown in Eq. p6p accounts for the 
collective light shifts and inhomogeneous broadening. In 
(38) this paper, we neglect the dipole shifts and frequency 
chirping by dropping the diagonal shift terms, and stress 
on the quantum correction that makes significance to the 
superradiance mechanism. The relevant part of the mas- 
ter equation now reads 
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We follow the derivation in [30|, |31| and summarize the 
main results. can be evaluated through 



2 



r = ^ 



r = ^ 



^2 



d3:r5-*(fo, t|f , A)5*-*(fo, i|f , A)p(i)^(r- , i| A) 

I 



(40) 



(41) 



We denote F = F^, and f = F^- for i ^ j- For above, 
the retarded function P'''^*^ the single-particle and two- 
particle source functions P^^)* and p(2)s ^j.^ gjygj^ pg_ 
spectively by 



pW«(fo,i|A) 
p(2)-(r'o,t|xi,X2,A) 



AAV 

;i2 



F(0 



(2+F)-zA 

2^(t)(i + r) 



(^+F)2 + A2 

2y(^)(i + r) 

(i+F)2 + A2 



(42) 
(43) 
(44) 



with M is the particle volume density, and — A accounts 
for the frequency of the Fourier component relative to 
that of the real field. Quantities A{t)^ ^(0; ^^'^ Y{t) are 
determined by the density matrix through Eqs. (|9ll0llip . 
In addition, The retarded Green's function is 



67reoc2 x 



(45) 



with q'a = uj/c and gp - 3^^^ 
integration for Eqs. (|40l4ip . we finally get 



&^pret(A)_ Through direct 



2C 



F = 



7(e 



F + 7/2 



V{t) 



2C'g 



2 4i'HCg) 



F + 7/2 



Y{t) (46) 



(47) 



with C(A) ^ Kg (r+,/2)/+A^ ^W, e~(A) ^ c.d/(2c) - 
(P^C(A), and /(C, ^) = [((C - l)e^ + cos gf + {ge'^ - 
sinp)^]/(C2_j_g2-j2 For large p and C, I{g) can be approx- 
imated by e2^/((:^2 _|_ g2y set A = for the resonant 
case. 

Note that C — CiCff) and the second term in Eq. I|46p is 
proportional toC2£i^e2^/(^2_|_^2^ ^ Q2g2^2C ^j^gj^ g'> C,- 

This indicates that the superradiant dynamics is, roughly 
speaking, characterized by the most relevant parameter 
Cg although the dependence on C and g individually can 
be evident, but less significant, through the exact form 
of Eqs. gSlHTD. 
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